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Abstract 

The low energy properties of the stationary state of a nonequilibrium Bose gas with the dynamical 
equilibrium between the escape and creation of particles are studied. The low energy spectrum 
of elementary excitations in such Bose gas is calculated at sufficiently low temperatures. Both in 
the nonequilibrium stationary state and in the thermodynamically equilibrium state the spectrum 
of the low energy excitations is formed by the long wavelength phase fluctuations of the complex 
Bose field. The spectrum of the phase fluctuations is found to have a diffusion behavior. Such 
behavior is in contrast to that in the thermodynamically equilibrium BEC systems in which the 
long wavelength phase fluctuations have the sound spectrum. The diffusion character of the phase 
fluctuations is due to noise generated by the escape and creation of particles. The diffusion character 
of the phase fluctuation spectrum results in the absence of coherency and, thus, in the absence of 
BEC and superfluidity as well. 



I. INTRODUCTION. 



The investigations of the coherent properties of a Bose gas of quasiparticles in the nonequi- 
librium systems attract large interest of experimentalists and theorists. The various efforts 
to realize the state with the Bose-Einstein condensate (BEC) in the nonequilibrium Bose gas 
of quasiparticles were carried out in different systems. The systems are, as follows, exciton 

no. 

gas in the double quantum wells in semiconductors PJ-pj], polariton systems, in particular, 
polariton systems in the two dimensional (2D) cavities located between two mirrors sj]- 15]. 
and a gas of excited magnons 16|- 19j. All these systems are the candidates to observe BEC 
at relatively large temperatures. Some of them are interesting both for the fundamental and 
for the applied science. In all these systems quasiparticles are created by external forces and 
live during a finite time. In what follows, quasiparticles are called as particles for brevity. 

The properties of the stationary state of a nonequilibrium Bose gas produced by the 
dynamic equilibrium between the escape and creation of particles due to external forces are 
considered in the present paper. In contrast to quantum atomic gases in traps in which BEC 



nas 
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3een observed at ultralow temperatures and where the atoms have the infinite lifetime 



21], all the systems considered here have a finite lifetime for Bose particles. Provided 



the lifetime is long compared with the thermalization time, the system can be observed 
within the time interval small compared with the particle lifetime and large compared with 
the thermalization time. In this case the external pumping is unnecessary. However, if the 
lifetime is small compared with the thermalization time, a continuous external pumping of 
particles becomes necessary to reach a stationary state of the system. Below we analyze this 
case within the framework of the nonideal Bose gas model augmented with the processes of 
the continuous escape and creation of particles. Moreover, the latter processes are supposed 
to be noncoherent in accordance with many experiments, for example, Q-Q], 10 1— 12]. In 



particular, for polaritons in microcavities, a pumping produces particles with the large 
energies and then the particles relax to the low energy states. We assume that the processes 
of the escape and creation of particles are governed by random sources. In the stationary 
state the escape and creation terms balance each other in average but the noise generated 
by these two processes acts on the system continuously. 

It is well known that the possibility of BEC and superfluidity is closely connected with 
the behavior of the low energy part of the excitation spectrum. The low energy spectrum of 
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elementary excitations is found here for the stationary state of a nonequilibrium Bose gas 
at sufficiently low temperatures. In the nonequilibrium Bose gas in the stationary state as 
well as in the thermodynamically equilibrium gas with BEC the renormalized low energy 
spectrum of elementary excitations is governed by the long wave phase fluctuations of the 
complex Bose field. As is well known, the low energy spectrum of phase fluctuations and, 
correspondingly, the low energy part of the excitation spectrum in the thermodynamically 



equilibrium Bose gas with BEC has a sound-like behaviour [22]- 24|. We show here that, 
in contradiction with the thermodynamically equilibrium Bose systems, the spectrum of 
long wave phase fluctuations for the stationary state in the nonequilibrium Bose gas with 
the noncoherent particle creation and escape has a diffusion character at sufficiently small 
temperatures. The diffusion behaviour of the spectrum is due to the noise generated by the 
escape and creation processes. This noise exists in spite of the balance in average between 
the escape and creation terms. The frequency behavior for the different components of the 
self-energy part of the one-particle Green function 29], 30], 25] determines the influence 
of the noise on the system. The relation between the kinetic and imaginary parts of the 
retarded or advanced components in the self-energy part gives an effective temperature of 
the stationary state. The important difference between thermodynamically equilibrium Bose 
systems and nonequilibrium systems in the stationary state results from nonzero value of 
the zero-frequency Fourier component in the kinetic term of the one-particle self-energy 
part of nonequilibrium systems. This circumstance does lead to the diffusion character of 
the low energy spectrum of the phase fluctuations in the nonequilibrium stationary state 
of a Bose gas. At the same time, the fluctuations of the Bose field module has the same 
character as in the thermodynamically equilibrium state of Bose gas. The diffusion character 
of the spectrum of phase fluctuations results in the absence of coherency and, thus, in 
the absence of BEC and superfluidity. In particular, for the systems with phase diffusion, 
the interference experiments should not demonstrate the pronounced interference picture 
in contrast to similar experiments in the quantum atomic Bose gases in traps 20], 21 ]. 
Unlike thermodynamically equilibrium Bose gas with BEC, the classical Bose field for the 
stationary state of nonequilibrium Bose gas obeys the Fokker-Planck equation which is more 
complicated than the Gross-Pitaevskii one. This is completely due to an existence of the 
noise generated by the particle escape and creation. So, if we try to describe the system by 
the Gross-Pitaevskii equation, the noise term should be augmented to the right-hand side 
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of the Gross-Pitaevskii equation. As a result, the equation goes over to the Langevin-like 
equation. 

At present, a majority of the systems used for the experimental BEC realization in 
nonequilibrium stationary state of quasiparticles have the quasi 2D geometrical configura- 
tions. For example, polaritons in microcavities are confined by two parallel mirror planes 
with a small spacing between them. In this case the bare spectrum of polaritons has a 
quadratic dispersion law in the region of sufficiently small momenta near the bottom of the 
spectrum. For the sufficiently small densities, such polariton systems can be described by 
;he nonlinear Schroedinger Hamiltonian with a noncoherent pumping and escape of particles 
8], {9]. The experiments on the BEC realization or laser effect in the polariton systems in 
microcavities are carried out actively at the present time 

0-Q. 

In the present paper the coherent properties and the low energy spectrum of excitations 
in the stationary state of nonequilibrium Bose gas are studied in the assumption that the 
system considered can be described by the nonlinear Schroedinger action with a noncoherent 
pumping and escape of particles. As will be seen below, in such systems the renormalized 
low energy spectrum has a diffusion character both in two and in three dimensions. To be 



closer to experiments 



10j - 12] in which the BEC of quasiparticles has been investigated, we 



consider the properties of a nonequilibrium Bose gas in the 2D geometrical configuration 
with a brief analysis for the typical properties of this dimensionality. 



II. LOW ENERGY SPECTRUM OF A NONEQUILIBRIUM BOSE GAS IN THE 
STATIONARY STATE. 

In the present paper we study the coherent properties of the stationary state in dilute 2D 
polariton gas in a microcavitiy under continuous noncoherent pumping and escape of parti- 
cles at small effective temperatures. These properties are determined by the particles with 
the momenta near the bottom of the polariton spectrum. The polaritons are concentrated 
in the region near the bottom of the spectrum after relaxation from the excited states of 
particles created by the noncoherent pumping. In the region of sufficiently small momenta 
the bare polariton spectrum has a quadratic dependence on the momentum for polariton gas 
in a microcavity confined by two plane mirrors separated by sufficiently small spacing |8[, 
9) . The renormalized spectrum of elementary excitations in the region of small momenta is 
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determined by an effective interparticle interaction. In this region of momenta the polariton 
system of sufficiently small density can be considered as a dilute 2D Bose gas with the ef- 
fective interparticle interaction independent of the momentum transfer. The Hamiltonian of 
such Bose gas with the omitted terms of the escape and pumping of particles can be taken 
in the form of the nonlinear Schroedinger Hamiltonian 

H = j d 2 r j> (-^ 2 -ty+\g (W) 2 } (1) 

where g > is the two-particle scattering amplitude independent of the momentum 
transfer, ijj are the creation and annihilation operators of Bose particles, m is the 
effective mass, and /i is the chemical potential. We put here H = 1. The particle density 
n =< > is supposed to be small. In the case of the low density Bose gas the chemical 
potential is fi = gn. The consideration of low energy excitations gives the possibility to 
take the bare spectrum of particles in the form p 2 /2m, or in the coordinate representation 
it can be written as ^— \^ 2 /2mj. The effect of the escape and creation processes will be 
included into the description of the system via introducing the relaxation terms into the 
Gross-Pitaevskii action treated within the framework of the Keldysh-Schwinger technique 



^{^t + v-^ 2 )^-\g (H) 

E (&f a + Ja*P)+7 a K a f< 



29], 30], 25] in the functional-integral formulation 31]. 

In the "tetragonal" representation of the Keldysh-Schwinger technique the action of the 
system S can be written in the terms of integrating over the double time contour with 
time reversion. The particle escape and creation processes can be included into action S as 
random sources of the particle escape and creation. The action can be written as 

(W 2 -l 

S=j>dtd 2 r{ x __^ (( - f ^-\ t ^- ~ f ^ ) (2) 

a=l,2 

where ip and complex conjugate ip are the Bose fields, /„, f a are random sources depend- 
ing on the time and space coordinates, index a = 1 corresponds to the particle escape and 
a = 2 corresponds to the particle creation by the external forces, and K a is the white noise 
correlator of random sources f a . The corresponding generating functional Z with the action 
S d2J) is written in the form of the functional integral over the fields ip and f a 

Z= j D^jD^ J] Df a Dj a exp {iS + i5S} (3) 



where the action SS depends on the infinitesimal sources ( as 



5S = j dtd 2 r (V< + C0) ( 4 ) 
The transition to the "triangular" representation can be performed by the substitution 



= -(^ + + ^_) 



(5) 



9 = i) + - ip. 



where ip + and ■0_ are the Bose fields at the upper and lower branches of the time con- 
tour. The component O is the "classical" component of the field and 9 is the component 
corresponding to quantum fluctuations. Note that the component is nonzero for the case 
of the fields coinciding at the different branches of the time countour ip + (t, "F") = ip_ (t, ~r^) 
, and the component 9 vanishes in this case. The generating functional for low energy par- 
ticles can be treated in the quasi classical approximation due to the large value of the low 
energy state occupation numbers. In this approximation the terms 99 can be neglected as 
compared with the terms GG. The transition to the fields Q, 9 and the integration over the 
random sources in the generating functional ([3]) with taking the inequality 99 « 00 
into account give the generating functional for the low energy fields. The corresponding 
action Sqci in the "triangular" representation [29|, |30| reads 



DQDQD9D9exp {iS QC i [9, 6] + iSS QC i [J,j]} 



(6) 



(7) 



Sqci [9, &] = j dtd 2 r <j +e ( idt + ^ _ -£_ + x f A ^\ Q 

- g ee (ee + ee) + i9v K (p) 9 

The frequencies and gradients of the fields describing low energy excitations and produc- 
ing the action Sqci [9, 0] are supposed to be small on the scale of \i = g < 00 >= gn, where 
fi is the chemical potential, ~$ = -iv, cD = idt are the momentum and frequency operators 
in the space-time representation. The frequency-dependent terms Vr(u), Ta(oj), T K (tu) 
are the imaginary parts of the corresponding self-energy part components of one-particle 
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Green function. The indices R, A, K denote the retarded, advanced and kinetic compo- 
nents, respectively. These terms can be expressed via the relaxation parts corresponding to 
the escape r* "*) (co) and creation r( m ) (co) of particles 



f R (co) = f g*> (co) - f £° (co) (8) 
f A (co) = T ( ° ut) (w) - f f ) (co) 
T K (co) = f %*> (co) + f K n) (co) 

For the stationary state, the relaxations Tr (co), Ta (co) vanish at zero frequency Tr (0) = 
Ta (0) = 0. In contrast, the term Tk (0) = (0) + Tjf 1 (0) is positive and nonzero since 



this term is a sum of two positive nonzero terms (0) and ?P (0) |29|. The kinetic 

component of the self-energy part Tk describes the noise generated by the particle escape 
and creation. For small frequencies co, or in the time representation for small time derivative 
co = idt, the relaxation terms can be expanded in co. The zero frequency components 
r^ (0) and (0) have zero value in the stationary state due to the dynamic balance of 
the particle escape and creation processes. For this reason, in the limit of small frequencies 
the relaxations Tr(co), Ya(co) can be expanded into a series in co. The first term of the 
expansion has the form 

f R (co) = -f A (co) = kco = xid t (9) 

where the constant x can be written as a ratio of (0) to some constant T* , which can 
be interpreted as an effective temperature of the stationary state of the system 

k=^t (10) 

ATP* v ' 

For small co, the dependence of Tk (co) on frequency co can be neglected as the zero 
frequency component of this term is nonzero Tk (0) = Tk ^ 0. The effective temperature 



T* is su 



pposcd to be much smaller than the Berezinskii-Kostcrlitz-Thouless temperature T c 



26], 27 1, 28 1 . Moreover, we suppose that T* << \i. The infinitesimal term 8Sqq} [J,j] is 



proportional to the infinitesimal sources J and j obtained from the infinitesimal sources ( 
by the transformation (J5J). In this case the sources £ + , £_ play a role of the fields tp + , ip_ 
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5S Q ci [J,j] = j dtd 2 r [(9J + J9) + (Oj + jO)] (11) 
Using notations ([9]) and ( fTUI) . we can rewrite the action Sqci [9, O] as 

Sqci[0,Q]= dtd 2 r{ 1 1 _* J _ > (12) 

where the functions F and F are 



^=1^+966-^16 (13) 

One can see that the functions F and F can be represented in the form of the functional 
derivatives 



F = dH ; F = dH (14) 

where d = 5/50 and d = 5/50 are the functional derivatives over the fields and 0, 
respectively. Thus, the functional Ho has the form of the energy 

H = J d 2 rO (j^-fi+ ^©©j © (15) 

The integral over the fields 9 and 9 in the generating functional Z (jBJ) has the Gauss form 
and can easily be calculated 

Z = jDeDOexp(-S[e)+i5S\j)) (16) 
The action S [0] can be written as 

S [O] = -!- I dtd 2 r {[i (1 + ix) d t O - F] [-i (1 - in) d t O -¥]} (17) 

L K J 

Note that, in the case IV — > due to the form of action ( TTTl) . the contribution to the 
functional integral f fl6|) is given by the configurations of the field obeying the standard 
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Gross-Pitaevskii equation. As is well known, the solution of this equation in the assumption 
of the small magnitudes of momenta compared with the inverse correlation length ^q 1 ~ 
^rrifi yields the sound spectrum of low energy excitations. 

In the case Tk 7^ the fluctuations of the Bose field should be taken into account in 
the integral for the generating functional ( fl6l) . Multiplying brackets in the action S [0] (117j) . 
we obtain 



S\ 



(18) 



= — Utd\ (1 + x2) ( ^ e) {dt ® ] + FF ^ 

v kJ y +i((d t e)F-F{d t e)) 

The complex fields and can be rewritten in the module-phase representation = 
pexp(i<p), = pexp( — icp). In this representation the corresponding terms of the action 
(fTBl take the form 



(a t e) (d t Q) = p 2 (d t <pf + (dtp) 



FF = fp*-±fp (^ 2 p-p(^) 2 ) 



+ 



where / is 



2m 



^p) 2 -2pN*p) (^Vf + 



( (V 7 ;)T • p-'(?vi 



+ 4((^p) (^))%4p((^p) (^))(^ 



(19) 
(20) 



f = gp 2 - n 



(21) 



and 



[(^e)F-F(^e)] 



(22) 



2fp 2 d t¥ > + 1 (<9 t p) [2 (^p) (^) + p^ 



pV 2 <p 



— (P<9i<p) 
m 



p-p 



The module p of the field can be represented as a sum of the average value and 
fluctuation of module p = p + 5p, where p is the bare value of the average module [28] 
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Po = a/ p/g, and dp is the module fluctuation. For small frequencies u « p, and momenta 
| p \« y/m/i the fluctuations 5p are small compared with_the average module 5p « p . 
The fluctuations of the phase in the 2D case can be large 
fluctuations 5p << po, the function / takes the form 



26(, |2Lj, 



. For small module 



(23) 



The expansion over the module fluctuations bp for the small values of the time and space 
derivatives (^^ 2 8p^ > (d t 5p) « pp leads to 

FF = W + 2p P0 Sp (I (^) 2 ) + pi (^^) 2 ) 2 + pi (^) 2 (24) 

i ( ($6) F — F (d t e)) = App Q 5p (d t cp) + 2p 2 (d t cp) (^~J ($<p) 2 (25) 



(d t Q) (d t Q) = p 2 (d t if) 2 + (dt6p) 



(26) 

The substitution of the above expressions into action ( IT8l) gives the generating functional 
and the corresponding action 



D5pD(p exp {iS [8p, ip] + i8S Sp ip } 



(27) 



S[6p,ip] = — I dt(f-r { 



+ 



(l + x 2 )pg(^ + 
4p 2 6p 2 + 4pp 5p(± 



• 2 f^^ 2 < 
^2m v 



(28) 



+App 5pd t¥ > + 2{p 2 d tV )^(^^ 
Integrating the generating functional Z (I2T|) over the module fluctuations Sp, we obtain 
the generating functional depending on the phase field (p alone 



Z = j Dip exp {iS [ip] + i5S^ [j] } (29) 

Here the effective action S [ip] for the small values of d t <p and j^-^ 2 tp compared with x/i 
if x < 1 or p, if x > 1 takes the form 



s M = ±- 1 dtd\ i x>pt (%>r + p 2 o ( ^VV ) } (30) 



r 



A' 




and 



6S V = J dtd 2 r [jpo exp (i<p) + jp exp (-i(p)] 
The action (130]) describes a diffusion. It can be rewritten as 

s[ v \ = ^f dtdh^f + (a,^V) 2 } (31) 



where Dn is the diffusion coefficient 



D ° = db (32) 



The action (l3Tj) gives the following expression for the phase correlator 



+ (d 2 

Near the pole this correlator can be rewritten as 



5 ( k - k')s(u-u') (33) 



< n jp-h. >= 7 — 4 r w x2p ° } \ ^ ("^ - 5 - ^ ( 34 ) 



The correlator describes the properties of low energy excitations. For the thermodynam- 
ically equilibrium Bose system, the phase correlator has the real pole of the sound type 
for small momenta and frequencies. In contrast for Tk ^ 0, the phase correlator ( jMJ) has 
the diffusion pole at the small momenta and frequencies. This form of the phase correlator 
means the absence of the phase coherence and, as a result, the lack of BEC and superfluidity 
in the stationary state of nonequilibrium Bose gas with the dynamic balance between the 
particle escape and creation processes. For example, in the case of the stationary state of 
nonequilibrium Bose gas the quantum interference experiments, like those in trapped quan- 



tum atomic Bose gases on the interference of two independent Bose condensates [20] , 21] , 
cannot give the pronounced interference picture. Note that the derivation of the low energy 
phase correlator fl33|) does not use the D=2 dimensionality of the system and can be used 
without any change for nonequilibrium D=3 Bose gas in the stationary state. 
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The phase correlator at the equal time moments can easily be obtained from Eq. (|33|) 

< Pt, t P-t,t >= ^k/^pD = ( 35 ) 

2_D /c 2 x Po fc 2 

The substitution of Eq. ffTOj) into Eq. ( 135|) gives the phase correlator in the form 

4T*m 1 , , 

< tPttf-tt >= —^1^ ^ 
Po k 2 

This form of the phase correlator is typical for the interacting Bose gas at low temperature. 
However, it results from the diffusion behavior of the excitation spectrum and not from the 
sound behavior of the spectrum. 

The transition to the space coordinates in the 2D case leads to the typical expression for 



the Bose field correlator at large distances and time moments 26j, 27] 



< @y,t@7V >= Po ex P \ ~~^~ ln (^) ) ( 37 ) 

where ^q 1 = a/ ragp\, s = y/\ r — r' | 2 +D \ (t — t') \ and ^ is the effective correlation 
length. 

It should be noted that the behavior of the system in the stationary state differs from 
the behavior of the equilibrium Bose gas in the Berezinskii-Kosterlitz-Thouless state. In 
the equilibrium Bose gas there is a local condensate with the coherent phase in the regions 
which size is much smaller than the average distance between vortices and larger than £o- 
The local condensate with the coherent phase in the stationary state of the nonequilibrium 
Bose system is absent. In this case the phase correlator has a diffusion character on all scales 
larger than £ - 



III. CONCLUSION. 



In the present paper the low energy properties of the stationary state of nonequilibrium 
Bose gas have been considered at sufficiently low temperatures T* << fi . In the stationary 
state the particle escape and creation compensate each other in average, i.e., the stationary 
state requires zero value of relaxations Tr (ou) and (u) for zero frequency u = 0, where 
^r,a (w) = T^"^ (u) — r^"| (u>). However, the quantum noise of the Bose field at zero 
frequency does not vanish and is characterized by IV = r^"^ (0) + (0), where the 
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kinetic component of the self-energy part IV at zero frequency is a sum of positive values 
p(o«t) ^ an( j to this noise the low energy spectrum of elementary excitations, 

which is governed by the phase fluctuations of the Bose field, has the diffusion character 
021]) • The diffusion character of the phase fluctuations means the absence of the phase 
coherency of the classical Bose field and, therefore, the absence of BEC in the stationary 
state of nonequilibrium Bose system. The diffusion character of the low energy spectrum of 
elementary excitations results from nonzero value of the zero frequency component of V k (w). 
Note that the behavior of the module of the Bose field in the stationary nonequilibrium 
state at small effective temperatures is the same as in the thermodynamically equilibrium 
state, i.e., module of the Bose field has homogeneous nonzero average value and the field 
fluctuations are small compared with the average value. Moreover, the momentum behavior 
of the correlator of the phase fluctuations in the nonequilibrium stationary state is analogous 
to that in the thermodynamically equilibrium state with the Bose condensate. 
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